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Abstract 

We develop an iterative algorithm for the genus expansion of the hermitian 
N x N one-matrix model (= the Penner model in an external field). By introducing 
moments of the external field, we prove that the genus g contribution to the m- 
loop correlator depends only on 3g — 2 + m lower moments (3g — 2 for the partition 
function) . We present the explicit results for the partition function and the one-loop 
correlator in genus one. We compare the correlators for the hermitian one-matrix 
model with those at zero momenta for c = 1 CFT and show an agreement of the 
one-loop correlators for genus zero. 
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1 Introduction 



Recently there has renewed interest in the hermitian one-matrix model which is defined 
by the partition function 



Z N [t.} = J DX e" (1.1) 



where the integration goes over N x N hermitian matrices. As is proven in Ref. 0, this 
model is equivalent as N — > oo to the following model in an external field: 

Z N [rj;a] = e-f^ 2 JDX e ^H^ X+ « log *) , (1.2) 

with 7] and {tk} being related by the Miwa transformation 

t k = ltir]- k - f5 k2 for k > 1, t = tr log rT 1 . (1.3) 

The partition function (|1.2p , in turn, is associated @ with an external field problem for 
the Penner model ||: 

Z N [A;a] = e-^^- 2 (detA)^) J DX e (A^^^) (L4) 

providing rj = A — aA~ x . The extra coefficients are introduced to provide 

Z aN [t] = Z N [r]; a] = Z N [A; a] (1.5) 

to any order of the genus expansion. 

A surprising property of the hermitian one-matrix model, which is advocated in 
Refs.Jl], [|, is that it reveals some features of c = 1 CFT interacting with 2D gravity. 
For the case of the Penner model, this property was discovered by Distler and Vafa |4j] 
in the double- scaling limit and has been extended to more general models by Tan ||, by 
Chaudhuri, Dykstra and Lykken and by Gilbert and Perry [0, in particular to the 
model (|1.4j) with A being proportional to a unit matrix. However, these features of c = 1 
have been observed in Refs.|l], 0] in genus zero and genus one, identifying a with the 
cosmological constant, without taking the double scaling limit. Therefore, it was conjec- 
tured that the hermitian one-matrix model has something to do with the continuum c = 1 
case similar to the Kontsevich model || which is associated with the continuum pure 2D 
gravity (c = 0) or to its generalizations || which are associated with c < 1. 

A direct way to verify this conjecture is to study correlation functions of the loop 
operators. For the hermitian one-matrix model with an arbitrary (not necessarily sym- 
metric) potential, the one-loop correlator has been known in genus zero for a long time 
flQ| , while the two- and three- loop correlators were obtained in Ref. [11]. The genus one 



contribution to the one-loop correlator was explicitly calculated in Ref.|]T2[ for a quartic 
symmetric potential and in Ref.pl] for an arbitrary symmetric potential. On the other 
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hand, much is known now about c = 1 correlators since pioneering works by Kostov [13 
and Boulatov [ I4"| . An incomplete list of references includes fl5|— [pi . 

In the present paper we develop an iterative algorithm for calculating the genus expan- 
sion of the hermitian N x N one-matrix model (= the Penner model in an external field). 
We introduce the moments, I p and J p , of the external field as well as the 'basis vectors', 
X*™) and ip( n \ which are determined by a recursion relation and diagonalize the iterative 
procedure. We prove that the genus g contribution to the partition function depends only 
on I p and J p with p < 3g — 2 while the genus g contribution to the m-loop correlator 
depends on I p and J p with p < 3g — 2 + m. We present the explicit calculation of the 
partition function and the one-loop correlator for an arbitrary potential in genus one. We 
compare the correlators for the hermitian one-matrix model with those at zero momenta 
for c = 1 CFT and show an agreement of the one-loop correlators (this correlator in c = 1 
is the only one which vanishes at non-zero momenta) in genus zero. We did not find such 
an agreement for the two-loop correlator. 



2 The iterative scheme 

We propose in this section a general iterative procedure for calculating higher genus 
contributions to the hermitian one-matrix model. We shall be solving, iteratively in 
1/iV 2 , the loop equation 

Ci 2m (A — uj) 
where 

oo r oo Q 

m = E^ fc , ^- = -EA-^ (2.2) 

fc=0 6V y A ) k=0 dt k 

and 

ww = W(X)^ m - W[y ' X) = WiX) W[X) (2 - 3) 

with the partition function Z[t] given by Eq. (|1 . 1|) . The one-loop correlator W(X) obeys 
the normalization condition 

XW(X) ^aN as A -> oo. (2.4) 



While we consider in this paper the loop equation (p7T|) , it should be noticed that a similar 



iterative procedure can be formulated for the equivalent Schwinger-Dyson equation 

where the partition function Z[rj_;a] is given by Eq.( |1.2| ) and {t^} are related to r)i - 
the eigenvalues of r\ — by Eq. (|1.3|) . The genus one solution of Ref. was obtained for a 
symmetric distribution of the eigenvalues by solving Eq.(|2] 
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Our goal is to solve Eq.( |2.1| ) iteratively order by order in 1/N 2 starting with the 
one-cut genus zero solution JTOl] 



W (X) = / ^ V(-) ^WT-c 



with b and c given by 



du V'(uj) r du uV'(lo) , T , n . 

- 0, / - . . == = 2aN. (2.7) 



'C*i 27TZ y/u 2 + bui + C iCi 27T2 ^Co- 2 + few + C 

In this way we will obtain the genus expansion of the (logarithm of the) partition function 

oo 

log Z[t] = Y J F g [t] (2-8) 

g=0 

and that of the associated multi-loop correlators. 
Following Ref.Q, we introduce the new variables 

I P = / 77- — r " 2 « + l)NS pQ , 

Cl 2m (CO 2 + ku + C f + 2 

J P = / $L YM ( 2 . 9 ) 

(c^ + ou; + c) 2 

Their relation to the moments of the external field rj, which is described by the density 
of eigenvalues p(x), can easily be established by substituting 



into Eq.(|J 



i/ p = /dx r-(| fo2 -| c )^o-5 P i, 

7V J (x 2 + 6x + c) p+ 2 

X J P = [dx ^ t-¥ 5 po- (2.H) 



^ 7 (x 2 + 6x + c) p+ 2 

The terms with <5 p o and 5 p i come from the second term on the r.h.s. of Eq. ( p.!0[) taking 
the residual at infinity. It is easy to see, expanding the integrands in Eq. ( |2.ll| ) in l/x 
that I p depends on with k > 2p while J p depends on with k > 2p + l. 

As has been proposed in Ref. depends at l < (7 < 00 only on J p and J p for 

p < P g — l (where P 9 = 3g — l as is proven below). This is in contrast to the t-dependence 
of F g which always depends on the whole set 

For an iterative solution of the loop equation, it is convenient to introduce the operator 
/CbyQ 

0>(A) = [ ^P^-Mu). (2.12) 

4 Inserting here Eq.( [2.l0| ), one finds the relation with the operator JC from Ref.|l[]: — jjK = K, + A, 
when applied to a function <£>(A) which decays at infinity faster than l/A (if < E'(A) ~ l/A as A — > 00, — l 
should be added to the r.h.s. which comes from the residual at infinity). 
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It is easy to calculate of how the operator /C acts on the set of functions 

$ (n) (A) = r (2.13) 

(A 2 + bX + c) n+ 2 

as well as on A^^^A). Let us start with n = 0. Comparing with Eq. (|2.6| ), one gets 
immediately 

/C$ (0) (A) = 2iy (A)$ (0) (A). (2.14) 

To calculate of how /C acts on & n \\) with n > 1, let us insert the following expansion of 
the denominator 

= (X + u + b)Y ~ — r — h -777; — fjr- for n>l. (2.15) 

A -co 1 (\2 + bX + c ) k (\ 2 + b\ + c) n X-cu ~ 1 ; 



Using Eq. (|2.9|) we get finally 

v ; id 2m (A - w) v ; ^ (A 2 + 6A + c) fc V ; V ; 

(2.16) 

Quite similarly, one calculates of how ft acts on A$^(A): 
/CA$ (n) (A) = -J n + A/C$ (n) (A) = 

n \ T — r T n ~ 1 T 

= E Z~l h uT +S m 7TV~ v> + 2Wb(A)A$<")(A). (2.17) 
fc=1 (A 2 + 6A + c) fc=1 (A 2 + bX + c) 

Some comments concerning these formulas are in order. The expansion is in integer 
powers of 1/ (A 2 + bX + c) due to Eq. ( [2.15| ). The sum over k is bounded from above by 
n by construction. For the symmetric potential V(A) = V(— A), which corresponds to a 
symmetric p(A) = p(—l) in Eq.( [2.1Q ), when 6 = and J p = 0, Eq. fl2.16|) recovers Eq.(5.13) 
of Ref.ll. 



Let us now turn to the iterative solution of Eq. (|2.1|) . Substituting the genus expansion 
for W(X) and W(X; A): 

00 00 
W(X) = Y,W g (X); W(X;X)=Y,W 9 (X;X), (2.18) 

and comparing the terms of the same order in 1/N 2 , one rewrites Eq. (|2.1|) as 

(lC-2W (X))W g (X)=Y,W g >(X)W g _ g/ (X)+W g „ 1 (X;X) (2.19) 

g'=i 

where g > 1. This equation allows us to calculate jy(A) order by order of genus expansion. 

As is proven below by induction, the r.h.s. of Eq. ( |2.19j) , which is known from previous 
orders of genus expansion, has a structure of the sum over n up to 

n g = 3g - 1 (2.20) 
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of the terms of the type A^/{\ 2 + bX + c) n and D^X/{X 2 + bX + c) n with the coefficients 
and being some functions of the moments I p and J p to be discussed latter. For 
this reason, it is convenient to introduce two sets of 'basis vectors' x (A) and ip( n \X) 
which obey 



/C-2^ (A))x (n) (A) 



(A 2 + bX + c) 



(/C-2^ (A))^ (n) (A) 



A 



(A 2 + bX + c) 



(2.21) 

This idea is analogous to the one suggested by Gross and Newman [25| for the unitary 
matrix model and for the hermitian model with a cubic potential. 

We can derive now recursion relations which defines x (A) and ip^^X) explicitly. 
Using Eqs. (l216D and ( |2~T7|) , we get 



X (n) (A) 



h - XJ 



A 



I$W(A) 



j n— 1 

XI X (fc) (A){il[/n+l-fc + &J n +l-fc] + Ji[cJ n+ i_jfc ~ J n -k]} 

k=l 
n-1 

^ {k \x){hJ n+l - k - jju+i-k] 

k=l 



A 
1 

A 



(2.22) 



and 

^W(A) 



A(ii + feJQ + cJi 
A 



$W(A) 



^ n— 1 

H X (fc) (A){cJl[/ n +l-fc + &«/n+l-fc] + (/l + bJi)[J n - k - Cj„, + i_fc]} 



fe=l 
n-1 



Y, ^ {k \X){ch J n+ i_ fc + (A + 6Ji)/ rt+ i_ fc } 



A 



fe=i 



where 



A = / 1 2 + 6/ 1 J 1 + cJ 1 2 . 



;2.23) 



(2.24) 



Eqs.fl2~2lD and ( |2~22D , ( [Og ) allow us to restore W 9 (A) provided the r.h.s. of Eq. (|219l) 
has the advertised form. The result reads 



W,(A) = £ [4 ra) x (n) (A) + Z^A^A) 



(2.25) 



n=l 



Notice that we do not add the terms with n = which are annihilated by the operator 
/C — 2140(A). These terms would contradict the boundary condition fl2.4p . The expression 
on the r.h.s. of Eq.( p.25|) is analytic everywhere in the complex A-plane except for the 
cut coinciding with the one of Wq(X). This is in accordance with the different iterative 
procedure advocated by Migdal and elaborated by David [pq*| . 

It remains to prove that the r.h.s. of Eq. (|2.19|) indeed has such a form. We see from 
Eq. fl2.25|) that the first term on the r.h.s. of Eq.( p,19[) has exactly this form order by order 
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in the genus expansion, while the second term needs a more careful analysis since we have 
to calculate 5W g -i(X) / 5V (A) . To do this we need some more formulas. 

From Eq.(|2.9|) one gets 



dl dl 

^ = -(p + l)I P+ i, ^ = -(p+|)^i (2-26) 

while Eq. ( |2.7p , which determines b and c, can be rewritten as 

I = -N, J = 0. (2.27) 

One needs as well the following formulas for calculating 5/5V(X): 

5b d 2 h-XJi 



5V(X) dX VA 2 + bX + c A 

5c d 2 (A + 6)/a + cJa 

5V{X) ~ dX^X 2 + bX+l A : 

with A given by Eq. (|2.24| ), and 

5I p d , , u ,\ , , x 5b , . , s r , , T . , , 5c 



(2.28) 



7^ A$(P) (A) - T777Tt(p + I) + + J P ] ~ T777Tt(p + D W 



^ - lr^HA)-TT#TT(P + l)Vi-77^(P+|)^ + i- (2-29) 



W(A) c)A v ' 5V(A)^ 2J 5V(X) 
These formulas^] can be obtained from Eqs.( ^9|) and ( |2.27| ) using Eq. (|2.26| ). 

It is easy to see that the result of acting of 5/5V(X) on W g -\ given by Eq.( [2.25|) has 



exactly the form discussed above. Moreover, we have presented an explicit algorithm for 
calculating W g (X), say by symbolic computer calculations with the only input parameter 
being g. The explicit results for genus one are presented in the next section. 

We perform now a power-counting analysis to express n g and P g , which are introduced 
above, via g. Let us note first that the highest moments Ip and Jp g emerge on the r.h.s. of 
Eq.( |2.25[) from the highest term 1/(A 2 + bX + c) n<3 on the r.h.s. of Eq. (|2.19| ) according to 
Eqs. fl2.22] ), ( |2.23| ) (they are associated with the k — 1 terms). Therefore, one gets P g = n g . 
To pass to the next order of the genus expansion, one analyze the structure of the r.h.s. 
of Eq.( [2.19j ). It is easy to estimate the highest power of each term W g ' (X)W > (A) for the 
solution (|2.25|) , which is known from the previous order, to be n g i + n g _ g > + 1 while the 
two-loop correlator gives the power P g + 3 by virtue of Eqs. Q2.28]) and ( 2.29 ). Therefore, 
one gets n g+ i = n g + 3 and finally 

P g = 3g - 1 (2.30) 



5 We have not calculated explicitly d/d\ in order that 8/5p(X), which is related to 8/5V(X) by 

S Id S 
8V{\) ~ ~N~dX6p(\) ' 

could be extracted easily as well. 
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since n\ = 2. Notice that both types of terms contribute to the maximal power since 
n g + n 9-a + 1 = 3ff - 1. 

According to Eqs.( |2.8| ), ( |2.3| ) and (|2.28| ), ( |2.29| ), the highest moments which contribute 
to F g are h g -2 and Jz g -2- An analogous result for the Kontsevich model has been obtained 
by Itzykson and Zuber P7|f| Applying Eqs. fl2.28]) , ( |2.29| ) m times to obtain the m-loop 



correlator, one finds highest moments to be I^ g -2+ m and J^ g -2+m- Thus, we have proven 
the following 

Theorem The genus g contribution to the partition function depends on 
I p and J p with p < 3g — 2. The genus g contribution to the m-loop correlator 
depends on I p and J p with p < 3g — 2 + m. 



3 The partition function and correlators in genus one 



We present in this section the explicit results for the genus one one-loop correlator Wi(A) 
and partition function Fx for the case of an arbitrary potential ^(A) which are obtained 
according to the algorithm of the previous section. 

In genus one there is no sum on the r.h.s. of Eq. ( [2.19| ) while |TT 



^o(A;A) 



b 2 -4c 



x2 ' 



16(A 2 + b\ + c 

Eq.( [2.21| ) yields immediately 

b 2 -4c 



2(A-^ 



A/i — |6(A + n) + c 



VA 2 + b\ + c + bfj, + c 



(3.1) 



^i(A) 



X (2) (A) 



16 



b 2 -4c ( $ (2) (A) 



16 



A 



[h-XJi]- 



$«(A) 



A 2 



[-J 2 h + l\l 2 + bI 2 J 2 - cJ\h + 2cJ 1 J 1 J 2 ] 



(3.2) 



where we have substituted the explicit form of X^ 2 \X) given by Eqs( [2.22D and ( p.23|) . For 
the symmetric case when 6 = and J p = 0, Eq.( p.2|) recovers Eq.(5.18) of Ref.Q. 

We are now in a position to integrate Eq. (|3.2|) (w.r.t. p) and find the functional F\ 
whose variation generates this W\{\). First, we integrate over A using the formulas 

6 It would be suggestive to relate these numbers to the (complex) dimension of the moduli space A4 g . m , 
which equals 3g — 3 + to. 
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1 / dx x + \b 



+ const, 



4 / (a; 2 + bx + cf' 2 (b 2 -4c)Vx 2 + bx + c 

+ 



1 / fix a; + \b 



4 -> (x 2 + bx + cf 2 " 3{b 2 -Ac){x 2 + bx + c) 3/2 

+ const. (3.3) 



\{x + \b) 



3{b 2 - 4c) Va; 2 + bx + c 

It is worth to note the appearance of the factor b 2 — 4c in the denominators. 

Using Eqs. (|2^9|) and ( gggD , we can now integrate W^.(A) explicitly. Namely, it is 
more or less clear that one can produce coefficients like in Eq. (|3.2|) by differentiating 
terms like log A (the terms with A 2 in the denominators originate from the variations 
( p.29| ), ( p.28| )). Moreover, due to the appearance of A( - fe 2 1 _ 4e ) term after the integration of 
lUi(A) over A, it is reasonable to assume the presence of log (4c — b 2 ) term in F\. After all 
these preliminaries, a straightforward, though lengthy, calculation gives the answer 

Fl = ~Y2 log(4c ~ fo2) ~ ki log A + const ' (3,4) 

which, in particular, reproduces Eq.(5.18) of Ref.jll for the reduced model (the symmetric 
potential) . 



4 Comparison with c = 1 correlators 

We compare in this section the explicit formulas for the loop correlators in the hermitian 
one-matrix model with those in c = I CFT. A problem immediately arises that the loop 
operator in c = 1 depends both on the length of the loop / and on momentum p which 
is associated with a matter field. Since in the hermitian one-matrix model there is no 
dependence on p, we should put p equal to some value (or sum up somehow over special 
values of p). The procedure is, however, unambiguous for the one- loop correlator which 
is non-vanishing only at p = due to the momentum conservation. Therefore, our idea is 
to compare first the one-loop correlators. 

The genus zero one- loop correlator, given by Eqs.( [2.6|) and ( |2.7|) , is greatly simplified 
after differentiating w.r.t. the cosmological constant a: 

l dW (X) _ 1 

N da VA 2 + b\ + c ' 

Let us put for simplicity 6 = keeping in mind that the general case can be restored by 
the shift A — > X + b/2; c — > c — b 2 /4. The equation (|2.7|) which determines c can be written 
in the form of the genus zero string equation 

If we restrict ourselves to the simplest case t 2 ^ 0, t 2m = for m > 2, Eq. (|4.2| ) gives 
c = —a/t2 and c > corresponds to the 'upside-down' potential which is familiar from 
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the quantum mechanical description of c=l correlators. Moreover, the r.h.s. of Eq. (|4.2|) 
is as well the known expression for the (diagonal) resolvent of the Schrodinger operator 
with the potential c(a) in genus zero. Therefore, we conclude that the genus zero one-loop 
correlators coincide. 

The genus zero two-loop correlator in c=l theory at zero momentum can easily be 
extracted, applying the Gegenbauer's addition theorem, from the results by Moore and 
Seiberg HQ: 

W (l uP = 0; l 2 ,p = 0) = K (^(/i + h)) (4-3) 

where K$ is the modified Bessel function of the third kind. 

It is easy to see that the Laplace transform of Eq.( |4.3j ) does not coincide with Eq.( pM| ). 
A possible way out (a pure speculative one) might be to sum up the general formula of 
Moore and Seiberg over the momenta which corresponds to the discrete special states 
of c = 1 CFT [111, 28, 29]. This problem deserves further investigation. 
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